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o : 

I Accurate QED calculations of the interelectron interaction corrections for the (ls2p)2^Pi, 

' (ls2p)2"^Pi two-electron configurations for ions with nuclear charge numbers 10 < Z < 92 are 

performed within the line profile approach. Total energies of these configurations are evaluated. 
, Employing the fully relativistic treatment based on the j-j coupling scheme these energy levels 

become quasi-degenerate in the region Z < 40. To treat such states within the framework of QED 
we utilize the line profile approach. The calculations are performed within the Coulomb gauge. 
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I. INTRODUCTION 



To provide accurate quantum-electrodynamical (QED) evaluations of energy levels for two- and three-electron 
"j^ I configurations of highly charged ions (HCI) become now an urgent problem in atomic physics. This can be explained 
^ by the growing number of experimental data and the necessity to use the energy levels for the evaluation of important 
Q . characteristics of HCI, such as e.g. transition probabilities and recombination cross sections. 
■ ^ ' In the past an approximate relativistic approach based on variational non-relativistic wave functions has been used 
for evaluating energy levels [1]. Numerous theoretical results for few-electron ions have been obtained within the 
' framework of fully relativistic many-body perturbation theory (RMBPT) and relativistic all-order many-body theory 
Oh' (AO) [2]. However, rigorous QED results, which allow for a consequent order-by-order improvement of the accuracy 
of the calculations become more and more desirable. 

The approximation of non-interacting electrons is commonly employed as a zeroth-order approximation in QED 
^ , perturbation theory for HCI in the absence of a quasi-degeneracy of levels. Accordingly, within the zeroth-order the 

■ energy of the few-electron configuration appears as the sum of the Dirac eigenvalues for the occupied one-electron 
' states. One-electron corrections (termed here as generalized Lamb shift) include QED radiative corrections, such as 
' the electron self energy (SE) and vacuum polarization (VP) as well as nuclear corrections, i.e., nuclear size (NS), 
' nuclear recoil (NR) and nuclear polarization (NP), respectively. Few years ago a status report of one-electron energy 
• corrections has been presented in detail in [3] . Since then considerable progress concerning the evaluation of higher- 

■ order self-energy corrections has been made [4] . 

, The dominant two-electron contribution is due to the interelectron interaction. Ab initio QED results for the 
^ ' first-order interelectron interaction in two-electron ions are known from Ref. [5] (see also [6]). The higher-order 
• ^ , corrections are much more intricate. Complete QED calculations of the second-order interelectron interaction have 
j>!Lj' been accomplished for the ground state and for non-degenerate low-lying excited states of He-like and Li-like ions 

^: [7-13]. 

The other important two-electron corrections are the screened self-energy and vacuum-polarization corrections, 
' which have been evaluated in [14-17] for non-degenerate two- and three-electron configurations in HCI. 

Various general bound-state QED approaches have been employed for the derivation of the energy corrections 
in HCI. The one most commonly used is the adiabatic S-matrix approach, developed by Gell-Mann and Low [18], 
^ generalized by Sucher [19] and first applied to bound-state QED corrections in Ref. [20] (see also [5]). The other one is 
■ - ' the Green's function approach, first developed in [21] and now applied frequently in a more sophisticated form of the 
two-time Green's function method [22-24]. Finally, the line profile approach (LPA) is utilized for the same purpose 
[25,26]. In our previous studies on the QED theory of interelectron interaction in HCI [11,13] this approach has been 
applied. 

In this paper we generalize the line profile approach to the case of quasi-degenerate electronic states in two- 
electron HCI. This problem arises, when a complete QED treatment including the relativistic j~j coupling scheme is 
applied to the fine-structure multiplets of systems with intermediate nuclear charge numbers Z . Some components 
of the multiplet with equal relativistic quantum numbers turn out to be close enough to each other (the (ls2p)2 ^Pi, 
{ls2p)2^Pi levels are the standard example). 
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Up to now the QED theory of the quasi-degenerate levels was considered only within the two-time Green's function 
method for the self-energy screening corrections (without any numerical applications) [27], for vacuum-polarization 
screening corrections [15] and within the covariant evolution-operator procedure [28] for the second-order intcrclcctron 
interaction. Numerical calculations for two-electron ions with Z = 10, 18 [28] are so far the only example where 
bound-state QED has been applied to the evaluation of the energy of quasi-degenerate levels taking into account the 
interelectron interaction up to second order. 

In this work we present an extension of the line profile approach, which is suitable for the evaluation of energies 
of any number of the non-degenerate or quasi- degenerate levels. The interelectron interaction corrections up to first- 
and second-order of QED perturbation theory are evaluated for the (ls2p)2 ^Pi and (ls2p)2 ^Pi levels in a wide range 
of Z values. 



The problem of the natural line profile in atomic physics was considered first in terms of quantum mechanics by 
Weisskopf and Wigner [29]. In terms of modern QED it was first formulated for one-electron atoms by Low [30]. In 
[30] the appearance of the Lorentz profile in the resonance approximation within the framework of QED was described 
and nonresonant corrections were estimated. Later the line profile QED theory was modified also for two-electron 
atoms [31] (see also [6,32]) and applied to the theory of overlapping resonances in two-electron HCI [33,34]. Another 
application was devoted to the theory of nonresonant corrections [35,36]. 

It was found in [25], that the LPA provides a convenient tool for calculating energy corrections. Moreover, it clearly 
determines the limit up to which the concept of the energy of the excited states has a physical meaning - that is the 
resonance approximation. The exact theorctic;al values for the energy of the excited states defined, by the poles in 
the Green's function, can be directly compared with measurable quantities only within the resonance approximation, 
where the line profile is described by the two parameters: energy E and width V. Beyond this approximation the 
evaluation of E and T should be replaced by the evaluation of the line profile for the particular process. Moreover, 
in the case of two-electron atoms the line profile approach was found to be very efficient for the evaluation of the 
reference state correction (reducible part of Feynman graphs) for two-electron atoms [37] . 



Consider the simplest process of photon scattering on a one-electron ion which is assumed to be in its ground state 
A (Fig. 1). Using the standard Feynman rules for bound-electron QED [6] yields the expression for the S-matrix 
element 



where V'a(^) = V'A(r)e is the wave function of the electron in the ground state, is the Dirac matrix together 
with the electron propagator 



A]j, ' {x) denotes the vector potential of the electromagnetic field (photon wave function). The notations Xu = {tu, r„) 
and Xd = {td,Vd) indicate "up" and "down" vertex coordinates, respectively. 

Insertion of the expressions for the electron propagator and the photon wave function yields 



II. LINE PROFILE APPROACH 



A. Line profile approach for one-electron ions 




(1) 




(2) 




n 



x[7^''4t'^)(rd)V'A(rd)]. 



(3) 
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Here oj = |k| and oj' = |k'| arc frequencies of the absorbed and emitted photons, respectively, k, k' arc the photon 
momenta and A, A' denote the photon polarizations. The summation over n is extended over the entire Dirac spectrum 
of electrons in the nuclear Coulomb field, e„ are the Dirac energy eigenvalues. Integrating over time variables {tu,td) 
and abbreviating the expressions in the square brackets by $y^(r„) and $^(rd), respectively, we can write 



27r -£„(!- lO) 

x$aM. (4) 

The function $^ (r) can be considered as a vertex function, which describes the absorption of a photon by an electron 
in its ground state. Below, we will formulate the resonance approximation, where we can define the energy and width 
of the level which have the general meaning independent of the features of the considered scattering process. Hence, 
the energy and the width will not depend on the function ^^(r) and thus we may consider the function as 
arbitrary. In particular, it can account for the interaction with the free electromagnetic field (radiative corrections). 
Let us introduce in Eq. (4) the matrix 



= i-e) J 



r Mr)^A{r) (5) 



and integrate over the frequency Employing the definition of the transition amplitude U via 

S = -2m 5{ijj - uj') U (6) 

we derive the expression for the amplitude 

We will consider the resonance case, when the frequency lo is close to the value w'^^^ = Sa — Sa + 0{oi), where a labels 
one of the exited states of an ion. In the resonance approximation we have to retain in Eq. (7) only the dominant 
term with n = a in the sum over n, i.e., 

C/(2) = ^laTaA ^ T* D-^T . (8) 

In order to simplify the application of the line profile approach to the many-electron ions we introduce the abbreviated 
notations 

T = TaA, (9) 
D=w- + , (10) 

= e„ . (11) 

Notice, that the function T describes the process of scattering. 

To obtain the Lorentz contour one has to insert the electron self-energy part in the internal electron line in Fig. 1. 
For simplicity we neglect the vacuum-polarization part. To the lowest order this leads to the graph depicted in Fig. 2 
and the corresponding expression for the scattering amplitude evaluated within the resonance approximation takes 
the form 



V^^\uj)D-^] T (12) 



with 

F«(a;) = e2 (ER(u; + e^)) . (13) 

Here T,ii{lu) is the rcnormalizcd electron sclf-encrgy operator. The upper index at the fimction V indicates the order 
of perturbation theory with respect to powers of the fine structure constant a for the graphs contributing to this 
function. Repeating these insertions in higher orders we can compose a geometric progression with the l-th term 
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Qi = u 



(2) 
Aa 



W - ea+SA. 



= T*D- 



T. 



(14) 



The resulting geometric progression is convergent for any lo except for values within the interval w G 
[So — £a — I^^"^' I) £a — £a |] close to the position of the resonance. Applying the formula for a convergent geo- 

metric progression one derives 



U 



Aa 



= ^T*D- 

1=0 



V^'\uj)D-^ 



I 

T 

T*T 



(15) 



Hence, the resonance is shifted into the complex plane and Eq. (15) is defined for all uj values on the real axis. Eq. (15) 
presents the analytic continuation of the expansion X)(^q Qi to the entire complex plane. 

Taking the square modulus of the amplitude (15), integrating over the directions of absorbed and emitted photons 
and summing over the polarizations we obtain the Lorentz profile for the absorption probability 



dW{uj) = 



a A 



27r (w + - Re {VW{io)} f + (Im {VW{u;)}y 



(16) 



Here dW (uj) is the probability for the absorption of a photon within the frequency interval uj, uj + du) and TaA is the 
partial width of the level a, associated with the transition a ^ A- 

Taking into account the correction depicted in Fig. 2 we improve the position of the resonance 



+1/(0) + Re [V^'\s, - Ea)} + 0{a^) ■ 



(17) 



Formula (16) defines the line profile of the process of scattering. Within the resonance approximation the line profile 
can be described by a Lorentz contour which is characterized by two parameters: the position of the resonance and 
the width. We define the energy shift for the state a as the shift of the resonance. The energy of the state a is 



E = o;-^ + = + Re {v'^'H^a - ^a)} + O(a') 



(18) 



and the width of the level as the width of the corresponding Lorentz contour at the position of the resonance 

r = -2 Im {Viio''^')} = -2 Im [v^^^s^ - + O(a') . (19) 

We note that the energy E of the level and its width T defined in the framework of the resonance approximation do 
not depend upon the function T (or $^) and, therefore, they do not depend upon the type of the scattering proc;c!ss. 
For example, we are free to consider not only the scattering of a photon but of some other particle as well, which 
couples to electrons. Going beyond the resonance approximation the line profile can no longer be described by a 
Lorentz contour and, consequently, the energy level can not be characterized only by two parameters E and F. In 
this case the evaluation of the energy levels should be replaced by the evaluation of the particular line profiles which, 
in general, depend upon the features of scattering process under consideration. 

The real part of the matrix element ( ER(eo) ] describes the lowest-order contribution to the Lamb shift and the 

imaginary part, which is finite and not subject to renormalization, represents the total radiative (single-quantum) 
width of the level a: 



AE: 



SE 



(SR(£a)) 



rSE 



(20) 



The other contribution to the lowest-order Lamb shift L^^ originates from the vacuum polarization. This correction 
does not contribute to the width Ta [6]. 

Studying the higher-order Lamb shift in one-electron atoms within the line profile approach, we have to account for 
the Feynman graph depicted in Fig. 3. For reason of simplicity, we will not consider the other second-order graphs. 

In the case rii = = a and n2 ^ a the graph in Fig. 4 can be viewed as a second-order self-energy insertion 
(loop-after-loop, irreducible part) in the graph Fig. 1 within the resonance approximation. We derive the following 
expression for the scattering amplitude 
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rr(6) _ rn 
^Aa — ^ Aa 



r(2) 



T*D- 



where 



_V / an \ /_? 



(21) 



(22) 



Note, that the singular term n = a is not included here by definition. This term was taken into account in the 
geometric progression described above and represents exactly the second term of this progression. Repeating the 
evaluations leading to Eq. (16) with 



(2) f V('\u;) + V(^\u;) 



= T*D- 



{V^^\ij)+V^^\(j))D-^ 



we obtain the improved resonance condition 

+ Re{ y(^) (CJ^°^) + (^rcs-, | _^ ^^^3^ _ _^res . 

Solving this equation for iv'^'^^ up to terms 0{a^) yields 

^res ^ _ ^y(O) ^ y(l) _ ^ y (2) _ ^^-j 



0. 



(23) 



(24) 



} + 0(a^). 



(25) 



The term V^'^\ea — £a) is the contribution of the irreducible part of the graph in Fig. 3. The derivative term 
corresponds to the reference state (n = a) correction. In particular, it coincides with the reference state correction 
that arises from the Fcynman graph in Fig. 3 after application of the adiabatic S-matrix method [38]. The other 
second-order electron self-energy (SESE) corrections are irreducible [38]. 



B. Line profile approach for many-electron ions (non-degenerate energy level) 

As in the onc-clcctron case wc will consider the process of photon scattering on an ion which is assumed to be in its 
lowest (ground) state. Investigating a non-degenerate energy level associated with a configuration containing at least 
one Is-electron (such as (1525)2^50, (ls2pi/2)2^fo; (1525)2^51, (ls)22si/2) (15)^2^1/2) etc) we can represent the wave 
function of the ground state via a proper combination of one-electron Dirac wave functions. A procedure based on 
this approach has been accomplished in [11,13]. However, for the investigation of quasi-degenerate levels or a doubly 
excited level the interelectron interac;tion corrections have to be taken into account in the wave function of the ground 
state. Here we restrict ourselves to two-electron ions. The generalization to N-electron ions will be presented at the 
end of this section. 

In the one-electron case we introduced the function $^ describing the process of scattering under consideration. To 
introduce such a function for the two-electron system we may consider first the simplest process of photon scattering 
on a two-electron ion disregarding the interelectron interaction corrections to the initial (ground) state. This process 
is depicted in Fig. 5, where the ground state is represented by two non-interacting electrons one of which absorbs (or 
emits) the photon. Accordingly, the ground-state wave function is given by 

^a{xi,X2) = ^det{tpis{xi)il>u{x2)} , (26) 

^a{xi,X2) = ^ dct{{l;uixi)'4'isix2)} , (27) 

where ipis{xi) = V'is(ri)e~*'^i"' *i , ipis{x2) are the Dirac one-electron functions with different projections of the total 
angular momentum. The bar over the one-electron functions indicates the Dirac conjugation. 
The S-matrix element corresponding to the graph in Fig. 5 can be written as 
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X j d'^Xui d'^Xu2 d'^Xdi d'^Xd^S^iVu^ - Yd^)5{tu^)5(td^) duJn^ ^ AiXm^Xu^)-^''"^ A*^^'^^^'\xui) 

j d'^Xui d^Xu2 d^Xdi d^xdj^ivm - rdJ6{tui)5{tdi) dun^ *^(a;„i,a;„j7''"2 A*(^fcJ"^'^(a;„2) 



+ 



"2 



(28) 



In order to employ the functions \E'^, \E'^ we introduced additional integrations d^x^ 2 d'^Xd-^ 2 '^'^(i"t(i,2 ~ 
^di,^) ^{iui,^) ^{^di,^)- The first and the second term in the square brackets represents graphs, where the photon 
is absorbed (emitted) by the first or by the second electron, respectively. Since the functions a{x\,X2), ^^{^11X2) 
are antisymmetric it would be sufficient to consider one of these terms only. 

As in the one-electron case we will look for the position of the resonance and employ the resonance approximation. 
It implies the neglect of the non-singular terms (evaluated at the resonance) in comparison with singTilar ones. The 
terms in the sum over ni, in Eq. (28) may contain a singularity at the position of the resonance only if they 
correspond to the positive-energy part of the Dirac spectrum. Accordingly, in Eq. (28) we can restrict to the terms 
with e„i > 0, £„2 > 0. 

Introducing the function $^(a;i, 2:2) as 

<i>^(a;i,a;2) -7^^4';^H-^i)*A(-^i,2;2)(5(ti -t2), (29) 
|.^(xi,X2) = ^A{xi,X2)r'A<!;'^^'\x^)5{U-t2), (30) 



we can write 



S^l^ = {-ief j rf*a;„id*a;„2 d^Xd-, d^Xd^ diOm dwn^ $^(a;„i,a;„2) 

i V'ni(r«i)V'ni(l'di) i •0n2(''tt2)V'«2(l'd2) 

2tt ^ UJni - £ni (1 - iO) 2tt ^ LOn^ - E^^ (1 - iO) 

x$^(a;di,a;dj. (31) 

Here we can employ the identity 

1 1 27r (5(tj„^-g„J 

(l^ni - £ni (1 - «0)) (w„2 - £n2 (1 " «0)) i (w„2 - £„2 (1 - iO)) 



-I \ (32) 

(-a;„, £„, -|-iO£„J(a;„2 - £712(1 - ^0)) 



which follows from the Sokhotsky formulas 



- = -m6{x)+V- , - = mb{x)+V-, - = —l,(x). 33 

a; -I- iO ^ ' X x-i^ ^ ' x a; iO -a; -|- iO i ^ ' ^ ' 

In view of the orthogonality of the Dirac functions and the asymmetry of the functions and the first term of 
Eq. (32) yields exactly Eq. (28). For £„j > the second term of (32) does not contribute when inserted in Eq. (31). 
As it was noticed above, for £„j < the second term does not develop any singularity at the position of the resonance 
and can be disregarded. 

Having performed the integration over the time variables {tui,tu2,tdi,td2) we arrive at 

S^2^ = {-ief{2-Kf j d^Yu^d^Vu2 d^Vd, d^Vd^ dum dw^^ $A(r„i,r„J 
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2tt ^ UJn, - e„i (1 - iO) 2tt ^ UJ^^ - Ens (1 - «0) 

X (5(i;^ +a; - w„i - a;„2 ) * A (rdi , rd2 ) , (34) 

where 

$A(ri,r2) =7''^4''^^(ri)*A(ri,r2), (35) 

$^(ri,r2) = *^(ri,r2)7''^4,'=''^'nri) , (36) 
Ea = Sis + Sis ■ (37) 

Formally the expression (34) is similar to the one in Eq. (4) which has been derived in the one-electron case. Taking 
into account interelectron interaction corrections to the ground state the function $^(ri,r2) and the energy Ea will 
become more complicated, in particular, the function $^(ri,r2) will depend on a;„j, Wn^- Nevertheless, the form of 
the expression (34) would remain unchanged. 

Below we will employ the resonance approximation defining the energy and width of the level such that they will not 
depend upon the features of the particular process of scattering. Since the function $^(ri, r2) carries all information 
about the process of scattering we can assume it to be arbitrary. 

Accordingly, it is convenient to introduce a graphical designation: a rectangle with a letter A inside (see Fig. 6). 
Lower and upper rectangles represent the functions $^(a;i,a;2) and $^(xi,a;2), respectively, which are defined as 

^a{xi,X2) = $^(ri, r2)e-**i(^^ +^^6{ti - ti) , (38) 
^a{xi,X2) = $^(ri, r2)e"i(^^ +^'H{t^ - t2) • (39) 

Here $^(ri,r2) denotes a complicated vertex function describing the scattering process under consideration, Ea is 
the energy of the groimd state, cj, uj' are the frequencies of the absorbed and emitted photons. 

We will look for the position of a resonance near w''®^ = E^'^'' — Ea + 0{a), where E^^^ =£„ + £{, is the energy of two 
non-interacting electrons. Applying the identity (32) to Eq. (34) one can see that the delta-function term is singular 
close to the resonance, while the other term remains regular (here we assume the function <I>^(ri, V2) to be arbitrary). 
The resonance approximation implies the neglect of the non-singular (at the resonance) terms in comparison with 
singular ones. Accordingly, within the framework of the resonance approximation the expression for the S-matrix 
becomes 

5^2) ^ (^-ief{2-Kf j d^Vu^d^Vu2 d^Td, d^u^ dwn^ dwn^ $A(r«i,r„2) 

X(5(w„i +U)n2 - EA-'^')5{'^m -Eni) 

"2 (''^2 j 



X ■ 



E 



2tt ^ W„2 - (1 - «0) 

ni,n2 

x5(i;^-|-a;- -a;„2)$A(rdi,rd2). (40) 
Integrating over , a;„2 in Eq. (40) and introducing the notation 

r„in2 A = (-e) y" d^nd^ra V^„,(ri)V;„,(r2)*^(ri,r2) , (41) 

we can express the corresponding amplitude for the scattering process in a form, similar to Eq. (7): 

f|(2) ^ ^■Arnw2^"i"2 A .^2) 

^ „t;r2'^-^"^-^"^+^^" 

Since we arc searching for the position of the resonance near uf^^ — E^^^ — Ea + 0(q:), we have to retain only terms 
in the denominator of (42) for which + ^^2 = E^°^ . We assume the energy level close to E^^^ = Sa + Sb to be 
non-degenerate and hence, within the resonance approximation the amplitude takes the form 

^(2) ^ TX T n,A ^ 2^+^-12. _ (43) 

£:f„^-E(°^ + EA 
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Here the summations run only over quantum numbers ni,n2, satisfying the condition +£712 = E^^h The matrices 
T and are given by 

(2^)riin2 ~-^ni«2Aj (44) 

{DUn, =i^- V^^°) + Ea (45) 

together with 

= . (46) 

As in the one-electron case T defines the type of scattering process under consideration. 

The interelectron interaction correction in first order is represented by the graph in Fig. 7. In order to evaluate 
this contribution which also shifts the position the resonance, one has to consider the graph in Fig. 8. In this paper 
we employ the Coulomb gauge together with the covariant metric. The photon propagator can be written as 

D'^:,M'^2) = ^J jni;'U\n\,n,)e-^^('^-'^^ , (47) 

where ri2 = |ri — r2| and 

Ik,Mri2)=^-^^, (48) 

The propagator D^^^^{xi,X2) corresponds to Coulomb photons, while D^^i^^{xi, X2) describes transverse (Breit) 
photons. The neglect of retardation implies the substitution -^^i^2(^'''i2) = -^^i/:i2(^'''i2)- We employ also the 
following notation 

/^'*(0)„,6'a6 = E / d'r^cfr2 I'^'^^Mri2) {^w {ri)^' M^i)) {^b' {r2h^' M^2)) • (50) 

The Lorentz indices /Xj should indicate that the Dirac matrices 7''* act on Dirac wave functions depending on variables 
Tj. The corresponding S-matrix element reads 

= (-«e)^ J d'^xi d'^X2 dQ d'^Xu^ d^Xu^ d^Xd^ d'^xa^ duui duJu^ dcodi du>d2 

^ ^a(^W1 ' ^U2 )^ ^ ^^(j^Ul ^112) 

2tt ^ ujui - £ui (1 - iO) 2tt ^ uuu2 - £«2 (1 - «0) 

i V'di(i-i)V'di(rdi) i V'd2(''2)V'd2(rrf2) 
2tt 4^ ujd, - ffdi (1 - iO) 2tt 4^ uJd2 - £d2 (1 - «0) 

x^/^.M2(I^^Ui2)e-"(*-*=) 

xe-''-^(^-+^^5{td^-td2)'i>Ai^d^,rd2), (51) 

while the summations over ui,U2,di,d2 run over the entire Dirac spectrum. Employing the identity (32) and retaining 
only terms which are singular near the resonance, i.e. keeping delta-function terms in (32) only, we obtain 

5^^^ = -27ri<5(a; - u')U^^^ 
= —2m5{w — u}') 



X 



. , Ea ~ ^ui — Su2 

UiU2did2 

XC /(I £di + £ui |)uiU2did2 

\ ^Tdrd2 A ■ (52) 

Ea +^ - ^di - ^d2 
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Within the resonance approximation we are left with terms satisfying the condition e^^ +£u2 = ^di +£^2 = E^^^ ■ 

The second-order interelectron interaction correction is represented by the graphs in Fig. 9. In order to apply the 
line profile approach to the contribution of the "box" graph of Fig. 9 we have to consider the graph depicted in 
Fig. 10 (a). The corresponding S-matrix element reads 

S^^ = {—ie)^ J d^xi d'^X2 d^x^ d'^X4 dfl dE d'^Xui d^Xu2 d^Xd^ d^Xd2 dcOui duJu2 d^d^ du)d2 
i V'«i(r«i)V'iii(i"3) i V'«2(''«2)V'm2(i"4) 



X 



2tt ^ uju, - £«i (1 - «0) 2tt ^ 'uJu2 - £«2 (1 - ^0) 



^ g-iWuj -t3)g-jw„2 (tu2 -t4)g-i0Jni (t3 1 ) g-iWnj (t4 -t2 )^(li3^(Ii4 

i V'ni(r3)V'ni(ri) i i/^n^ (r4)-0„2 (r2) 



X ■ 



,,1 _ e-„ (^ — Ott 2^ , 



'^"1 - (1 - «0) 27r ^ - e^^ (1 - iO) 

X e~*'^'^l "'''l e~^'^''2 -*<i2 )'yA'l /yM2 

^ V'di(ri)V^di(rdi) i V'd^(r2)Vid2(i'd2) 
27r 4^ Wdi - - iO) 27r ^ ujd2 - £d2(l - «0) 



X 



x^J^,^,(|S|,ri2)e--(*-*^)^/^3M,(|l^|,r34)e-"(*3-t4) 

X e-^*<<i +-) - td2 {rd, , rd2 ) , (53) 

Employing the identity (32) and retaining only the terms in the summation over 1*1,1*2 and ^1,^2 which become 
singular close to the resonance, we obtain the following expression for the S-matrix element corresponding to the 
graph in Fig. 9 (a) 

Sf = -2nid{L0 - J)U^^^ 
= —2m5{LO — Lo') 



X 



, , Ea — ^ui — £«2 

«i«2di<i2 



Xg4 ^ f JO -^(|^|)'UiM2nin2-^(l ^ ^di + | )raira2c!i(i2 



2^ i (-1^ -f- £„i - £„i (1 - iO))(£;^ -t-w -t- n - £„i - £„, (1 - «o)) 

Tdid2A- (54) 



nin2 
1 



Ea +w - £<ii - £<i2 



Again within the resonance approximation only terms satisfying the condition £„j -|- £„2 = Sd-^ -\- 6^2 = E^^^ will be 
kept in the summations over ui,U2,di,d2- 

Let us consider separately the reference states terms, i.e., for which £„j -|-£7i2 = E^°^ holds. Inserting a similar 
identity for the energy denominators 

1 



(-0 + £„, - £„, (1 - iO)){EA +UJ + n-£m- Sn2 (1 " '^O)) 
_ 277 S{n- S^^+SnJ 
i (EA+t^-eni -£^2) 

_l zl (55) 

(O - £„j -I- £„j -h £„i iO)(i;^ -Fw O - £„i - £„, (1 - iO)) 

into Eq. (54) , one can verify that the term with the delta- function coincides with the second element of the geometric 
progression for the graph in Fig. 7. Hence, while generating the geometric progression this term will refer to the 
second element of the progression. 

In order to evaluate rigorously the position of the resonance up to second order in a we have to consider all 
corrections of first and second order simultaneously. Up to first-order of perturbation theory we have to account for 
self-energy (SE) and vacuum-polarization (VP) corrections as well as for the exchange of one Coulomb or one Breit 
photon. In second order we have to account for all one- and two-electron Feynman graphs of second order including 
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radiative corrections, screening of the self energy and vacuum polarization and two-photon exchange graphs. However, 
the evaluation of the radiative corrections is not the goal of the present work. Below we will present the derivation 
of a formula for the one-photon exchange contribution in Fig. 7, the two-photon exchange "box" graph on Fig. 9 (a) 
and the three-photon exchange "box" graph Fig. 11. 
The scattering amplitude can be written as 

U^^l = T+D-^ [(1^(1) + + y(3))£,-il T, (56) 
where V^^^ corresponds to the one-photon exchange graph Fig. 7 (see Eq. (52)) 

y« I^{\-£a + ea'\)a'b'ab. (57) 

g=c,t 

In contrast to the one-electron radiative corrections (see formula (13)), this one-photon exchange correction does not 

depend on lo. 

Taking into account the "box" graph in Fig. 9 (a) we obtain (see formulas (54), (55)) 

oo 

gg'=c,t £„i -I- e„2 jtea + £i,_oo 

1 



{-n + Ea' - e„i (1 - iO)){EA +LO + fl-ea'- (1 - «0)) 

oo 

(58) 



gg'=C,t £ni + Sn2 =^a "I" Sb_ 

-1 



(O - -he„i ?0)(£'^-Fw + n-ea' -en2(l - iO)) ' 

Again the summations over g, g' run over scalar (Coulomb) and transverse (Breit) photons, respectively. The second 
term in Eq. (58) represents the remainder after subtracting off the reference state singularity. This subtraction was 
done at the stage of generating the geometric progression with the one-photon exchange insertion. In particular, the 
identity (55) was employed for deriving the expression for the reference state contribution (the reducible part) of 
the "box" graph, where the delta-function term coincides with the second element of geometric progression for the 
one-photon exchange graph. 

In the case of one-electron ions the reference state term (see Eq. (22)) appeared only via the derivative term in 
Eq. (25). Considering photon exchange in two-electron ions the situation is different. In this case a reference state 
contribution appears directly in V^^\uj) (see Eq. (58)) while a derivative term does not arise since V^-^^ does not 
depend on ui. Nevertheless, if the sets {a, b} and {a', b'} are equivalent, the term corresponding to the reference states 
can be expressed as a derivative. At the point of the resonance we can set u) = — Ea + E'^^^ in Eq. (58). Accordingly, 
both factors in the denominator will be identical. Utilizing the formula 



{x + iO)2 dx [x + iO) 

and integrating by parts we can shift the derivative to /(Q). 

Let us now turn to the three-photon exchange correction neglecting retardation effects, crossed-photon graphs and 
the contribution of the negative energy part of the Dirac spectrum. Within this approximation the photon propagator 
does not depend on the frequency 51, which yields 

V / / y / y a'b'nzm n3rnnin2 n\n2ab 

gg'g"— c,t nin2U3n4 

(60) 



(£„3 -I- - Ea -w)(e„i + £„2 - Ea -w) ' 



where the prime at the summation symbol indicates that the reference states (£„i -|- £„2 = £a + £6, £713 + £714 = £« + £&) 
are omitted. 



10 



Taking together the contributions of V^^\ V^"^^ and V^^\ we can generate a geometric progression, where the Z-th 
term reads 



T. 



(61) 



Performing similar steps as in the one-electron case we sum up this progression and derive a condition for the position 
of the resonance 



+ Re{ V^^\io'^^) + V^^\ij'^^) + } - £^ = . 

The energy and the width of the level will be equal to (see Eqs. (18), (19)) 

-2 Im{ V^^^uj"'') + V^'^\uj"'') + V^^\ij"'^)} + 0{a^) . 



(62) 



(63) 
(64) 



It is important to emphasize that Eq. (64) has the meaning of the correction to the width of the level only if the 
full set of Feynman graphs of a given order is under consideration. Indeed, the contribution of the graph in Fig. 7 
cancels completely with a part of the contribution of the self-energy correction, while the vacuum polarization gives 
zero contribution to the width. Such a cancellation is an immediate consequence of the Pauli principle according to 
which transitions of electrons into occupied states are prohibited [39]. 

We note that the two-electron graphs to the first order in V'(^) does not depend on u>. Hence, the solution of Eq. (62) 
together with F^i), V^'^\ V^^^ given by Eqs. (57), (58), (60) yields 



Ea + Re{ y(i)(^W - Ea) + V^^\e^''^ - Ea) 



+y(3)(^(0) 
+ 0(a^). 



■Ea)+V^'\e^'^-Ea) 



dv^^\w) 



^ = E^^)-Ea 



} 



(65) 



The term V^^\ea — Sa) represents the contribution of one-photon exchange graph Fig. 7, while the term V^'^''{£a — £a) 
accounts for the contribution of the two-photon exchange graphs in Fig. 9. In particular, this term includes the 
contribution of the reference states occurring in these graphs. The third-order term V'^'^\ea — e^) contributes to the 
throe-photon exchange graphs in Fig. 11. It does not contain the contribution of the reference states because we 
disregard the retardation effects considering it within the framework of relativistic many-body perturbation theory 
(RMBPT). The derivative term in Eq. (65) as well as that term in Eq. (25) docs not correspond to a certain Feynman 
graph. Similar to Eq. (25) it can be related to the contribution of the reducible part (reference states) of the graph 
in Fig. 11. 

Let us mention that if we would take into account in addition the radiative corrections and as well as screening 
effects, we would obtain a corresponding correction 1^'^' (uj) containing the contribution of the electron self energy 
(13), the one-photon exchange (57) and the vacuum polarization. For V^'^\uj) we would similarly obtain the sum of 
Eqs. (22), (58) and in addition all the missing radiative effects of second order together with the screened self-energy 
and vacuum-polarization corrections. Accordingly, instead of Eq. (65) we would have been left with 

= -Ea +V^°^ + Re { V^^Ke"^^^ - Ea) + V^^\e^''^ - Ea) 
"5y(i)(w) 



+y(l)(£;(0)_^^) 

+\v('\e^'^-Ea) 
+V^^){E^o)-Ea) 
+V^'\e^°^ - Ea) 

+ y(2)(£;(0)_^^) 



du! 



a2y(i)(w) 



duj 
■5y(i)(w) 



UJ = ijCO) - Ea 
u) = £'(0) - Ea 
UJ = - Ea 



duj 



w = ijCo) - Ea 



+ 0{a^) . 



(66) 



Formulating the line profile approach for A''-electron ions, it might be convenient to introduce the function 



11 



N 



j=2 



(67) 



which should be depicted graphically by a rectangle with a letter A inside and with N out-going electron lines. Here 
• • • ,rjv) describes the TV-electron ions in the lowest (ground) state A together with the absorbed photon. 
Accordingly, formula (32) generalizes to 



JV 

n 



1 



\ (w„,. - £„^(1 - iO)) 



This identity can be written as 



N-l 

n 



— ()(CJ„ -£nj) + 7 — , -r,^ ^ 



1 



{oJnti - - «0)) 



(68) 



JV-l 



N 

n 



n f<5(a;„,-£„J 

\ {uirij - En,- (1 - ^0)) (u;„„ - £„„ (1 - iO)) 



f{uj,s). 



(69) 



Then the terms in Eqs. (34), (51) and (53), which correspond to the function f{w,e) will not contain the singularities 
and will be omitted within the framework of the resonance approximation. Eqs. (65) and (66) will remain unchanged, 
however, now V will contain additional contributions of three- and up to A^-electron graphs. In particular, for 
three-electron ions the functions ^^^^(0;) and V'^^^u) will also account for contributions of three-electron graphs (see 
[11])- 



C. Line profile approach for many-electron ions (quasi-degenerate energy levels) 

We now turn to the application of the line profile approach to quasi-degenerate levels. Without loss of generality, we 
can restrict ourselves to two mixing configurations. We will search for the positions of the resonances corresponding 
to these configurations and will construct basic wave functions ^'i and ^'2 within the j-j coupling scheme. The 
energies corresponding to these wave functions are denoted by Ei'\ and they are supposed to be close to the 
exact energies of the electron configurations under consideration. Employing the line profile approach we will consider 
a scattering of a photon on a two-clcctron ion in its ground state A- The positions of resonances may be found near 
the values wj*'*' = E^'' — Ea + 0(a) and lo'^''^ = £2*^ — Ea + 0{a), respectively. Within the resonance approximation 
we will have to retain two terms in the sum (42) corresponding to the basic functions ^'i and ^'2- The scattering 
amplitude may be written as 

UAa = T+D-^ [\VD-^] T, (70) 
where £> is a matrix 2x2, defined on the functions ^i, ^2' 

D =lu + Ea-V^''K (71) 
1/(0) =hi + h2, (72) 
AV = V- V^°^ = y(^) + 1/(2) + y(3) + . . . . (73) 

Here hi, /12 are the one-electron Dirac Hamiltonians acting on the one-electron Dirac wave functions depending on 
ri or r2, respectively. Since the functions and ^2 are orthogonal the matrix D is diagonal. Accordingly, we now 
have to compose a geometric matrix progression with the Z-th term 

Qi = T+D-'^ [AVD-'^] ' T (74) 
and sum it up employing the formula for a convergent geometric progression. The expression for the amplitude reads 

Ua = T+ [D - AVr^ T = T+ —^——T = T+ -T . (75) 
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Introducing the function $ = ($i, $2) by means of the relation $ = B^, where the matrix B is assumed to diagonalize 
the matrix V = V''°^ + AV, i.e. V^'""^ = B+VB. The expression for the amplitude can now be written in the following 
form 

~ -^A $1 , , I rR+T/Ri -^^i A -^A 4>2 ,,,17, Tr+tTri '^'''^ ^ 

= , .diag / s "^^i A + T'^is — , ,diag f .'^^2 A- (76) 

Taking a square modulus of the amplitude (76) and integrating over the directions of the absorbed and emitted 
photons yields a line profile for the probability of photon absorption. The positions of the resonances are determined 
by the equations 



ujr + Ea- Re{0' (^r)} = , (77) 
+ Ea- mV^lZ = • (78) 

Hence, the energy of the configurations are: 

E^,=Re{V^%\icjr)}, (79) 
E^,=Re{V^%\iLor)}. (80) 

Assuming that the energies of the configurations are close to each other, we can expand Eqs. (77), (78) into a Taylor 

series around the values ujI'^^ = — Ea + Ei^^ and lo^'^^ = — Ea + i?2°^ ' respectively. As in the case of non-degenerate 
levels this can be achieved up to any desired accuracy. 

Note, that employing the resonance approximation in case of non-degenerate level we have to retain in a correspond- 
ing sum (42) certain many-electron functions composed within the j-j coupling scheme. Indeed, after diagonalization 
of the matrix V all other combinations of one-electron functions will yield zero in view of the antisymmetry of the 
wave function of the ground state and the symmetry of the matrix V. Hence, having constructed a many-electron 
function in the j-j coupling scheme, Eq. (56) becomes a scalar one. 

The line profile approach outlined above can be easily employed for an arbitrary number of degenerate levels. The 
generalization of the method to AT-electron ions was described at the end of the previous section. 



III. EVALUATION OF THE ENERGY LEVELS OF QUASI-DEGENERATE TWO-ELECTRON 

CONFIGURATIONS 

We will evaluate the interelectron interaction correction for the two-electron configurations (ls2p)2^Pi and 

(ls2p)2^Pi. Employing the relativistic j j coupling scheme these energy levels become quasi-degenerate in the region 
Z < 40. To treat such states within the framework of QED we will apply the line profile approach. Within the j-j 
coupling scheme the wave function of a two-electron configuration can be represented as 

*JMju2iii2(ri,r2) = Af J2 CjM(wim2) 

mi, 7712 

X bl^jihrmiri)tpj2hm2{r2) - V'ji;imi(r2)V'j2;2m2(ri)] , (81) 

where the normalization constant is A'' = 1/2 for equivalent electrons and A'' = l/-\/2 for non-equivalent electrons, 
respectively. Cj^j (TOim2) is a Clebsch-Gordan coefficient. The one-electron Dirac functions V'j;m(r) are characterized 
by the standard set of one-electron quantum numbers: total angular momentum j, its projection m and the orbital 
angular momentum that fixes the parity of the state. For the two-electron wave function the relevant quantum 
numbers arc the total angular momentum J and its projection M. 

Following the procedure described in Section II C we will construct the matrix V (73) on the functions (81) 

*J=l,M=0,ii = l/2,j2 = l/2,ii=0,i2 = l = (l''^2pi/2) , (82) 
^'j=l,M=0,ji = l/2j2=3/2,ii=0,i2 = l = (l'5'2p3/2) (83) 

and examine the positions of the resonances close to loI'^^ = — Ea + ^is + £2pi/2 + 0{a) and Wj**^ = 
- Ea + Sis + £2p3/2 +0{a), respectively. 
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As it has been elaborated in Section II the operator V, in general, depends on lo. The position of the resonance 
can be derived via Taylor expansion at the approximate positions of the resonances LjjY'^ = — Ea + £is + ^2pi/2 ^^'^ 
'^'2^ = — Ea + £i« + £^2p3/2 E<1- (65)). For the practical calculations it is convenient to expand some matrix elements 
of V at the point ujI°^ and others at the point UJ2^^ keeping only terms 0{ct^) in both expansions. The resulting 
inaccuracy can be referred to corrections 0(q;^) [24], because at low Z values the energy difference £2^1/2 ~^2p3/2 
becomes small, while at large-Z values the degeneracy of the levels 2^Pi, 2^Pi is nearly negligible. 

The interelectron interaction correction is represented by the set of graphs (Figs. 7 and 9) which is symmetric under 
interchange of the upper and the lower indices and relabeling of the electrons in the graphs. Accordingly, the operator 
V is given by a symmetric (and in general complex) matrix. However, as a consequence of performing the Taylor 
expansion of the matrix elements of V and neglecting third- and higher-orders terms it can lead to a non-symmetrical 
matrix. To prevent this asymmetry arising due to purely technical reasons one may symmetrize the matrix V by 
hand. 

Hence, the matrix elements of the frequency-dependent operator V evaluated at the resonances can be written as 

{{ls2p,/2)\V{cj'n\{^s2p,^2)) = ((l.s2pi/2)|F|(ls2pi/2)) , (84) 

{{ls2ps/2)\V{u''^')\{ls2py2)) = ((ls2p3/2)|F|(ls2p3/2)) , (85) 

((ls2pi/2)|V^(a;--^^)|(ls2p3/2)) = ^ [((ls2pi/2)|i^|(ls2p3/2)) 

+ ((ls2p3/2)|i^|(ls2pi/2))] , (86) 
((ls2p3/2)|y(a;'--)|(ls2pi/2)) = ((ls2pi/2)|F(a;'--)|(ls2p3/2)) . (87) 

The operator F is defined via its action on the set of the one-electron Dirac functions {ab}, which in our case consists 
of {ab} = {ls2pi/2}, {ls2p3/2}. To zeroth-order perturbation theory the operator F reads (see Eq. (72)): 

F!,?lab=^aSa',a+ebSl,',b. (88) 

Being interested in ionization energies it is more convenient to introduced a shifted F with the zeroth-order matrix 
element 

Fa"b'ab = ^o. Sa',a + ^b h' ,b - ^Is -m , (89) 

where the electron rest energy m (in relativistic unites) and Is-electron energy are subtracted. In first-order per- 
turbation theory the interelectron interaction represented by the graph in Fig. 7 can be described by the matrix 
element 

Fi'lab = ^^Ii\^a'-ea\)a'b'ab. (90) 

Since graph in Fig. 7 is irreducible F^j^i^j^ coincides with the expression (Eq. (57)) for non-degenerate levels. In 
second-order perturbation theory we have to account for the two-photon exchange corrections depicted in Fig. 9: 



gg' nin2 

oo 

IHM)a'b'n,nJ^ (|0 - + 



I {e!°^ - Ek°L ) - + - +iO £„2 ) 



+ 1. /^^ /nN).W.../Ml»-^»+^.-|)....a 
^^J {E^°^-Ek°Lm-Sr,,+Sa'+iOen,) 



p(2)(box,red) 
a'b'ab 



gg' nin2 



r<0) 



(0) 



OO 

2^ / 



I^{\^\)a'b'nin2l^ (1^ ~ ^a' + Dnmzab 
(0-£„2+46^-£a' +iO£„2)2 
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27r 



dfl 



/S (|f2 - Ea + Sa' \) n\niba 



+ 



p(0) _ p(0) 
Jlinin2 



(92) 



p(2) (cross, irr) 4 J /l X ^ 

gg' nin2 ^ 

OC- J 

1^ ji^^-^qI^ Ea' ~\~Ea |)nia'6n2 



X— dfl 



277 7 (Ens -Sni +£6 -£a')(^ - +£a+«0e„2) 

— oo 

+ (1 - <5o,(£„2 -e„i -£(, + £„')) 

oo , 

(e„2 - £„i - £6 + £a' ) (f^ - £n2 + -E^6^ " ^a' +«0 £„2 ) 



2n 



(93) 



p(2)(cross,red) _ „4 c 

gg' nini 

oo , 

■ 27r y (O - £„2 +£a +^0£„J 



In Eqs. (91) and (92) the notations E^jj = £a £'I''b' = ^g' +£b'j £^n°ri2 = £rii + £ri2 ^ii'C introduced. Index g runs 

over 0, t (scalar and transverse photons). The Kronecker symbols ensure that terms with potentially zero denominators 
will be omitted in the summation over n\,n2- 

Note , that the Eqs. (91), (93), (94) for the irreducible parts coincide generically with Eq. (58) for non-dcgcncrate 
levels (see also [13]). However, for the reducible part of the "box" graph additional terms originating from the 
geometric progression for the one-photon exchange graph (non-diagonal matrix elements of the second term of the 
progression) occur. It is easy to make sure that the contribution of reference states {eIi^u^ = E^^J , Enin2 = E^^ili) to 
the exchange of two Coulomb photons (or Breit photons with neglect of retardation) is absent. 



IV. NUMERICAL RESULTS AND THEIR ANALYSIS 



The results of the numerical calculations are presented in Tables I, II. To account for nuclear size corrections we 
solved the Dirac equation with the Coulomb potential generated by a nuclear charge density described by a Fermi 
distribution. The parameters of the Fermi distribution are taken from Ref. [13]. 

In Table I we present a detailed analysis of our results obtained for the photon-exchange contribution. The value 
is the binding energy of the 2p state according to Eq. (89), the value V^^^ corresponds to the one-photon exchange 
contribution Eq. (90) and V^"^^ represents the two-photon exchange contributions given by Eqs. (91 94). We note, 
that, in general, the matrix V has complex elements and both their real and imaginary parts contribute to the energy 
eigenvalues, i.e., the real part of the diagonalized matrix V. In our calculation the imaginary part of V'^^^ is taken into 
account, while the imaginary part of the two-photon exchange contribution [V^'^^) is neglected. The values Eyi{2^Pi) 
and £'x(2'^Pi) denote the photon-exchange contribution to the energies of the corresponding electron configurations 
(neglecting the radiative corrections). For Z = 10, 18 we present also the values for the difference between the energies 
of the levels under consideration reported in Ref. [28]. 

In order to analyse the influence of the quasi-degeneracy on QED effects in more detail we compile the corresponding 
energy shifts of the levels due to the photon-exchange contribution calculated within various approximations. The 
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differences between the energies of the levels ealeulated without the approximations and the energies calculated within 
the framework of the approximations are presented in Table I [AEx{2^'^Pi)- Appr. 1-5]. 

Approximation 1: We omit the non-diagonal elements of the matrix V = V^^^ + V^^^ + V^'^\ Consequently, effects 
of the quasi-degeneracy are totally neglected. 

Approximation 2: We omit the non-diagonal elements only in the matrix V^'^\ As stated above the expression for 
the one-photon exchange correction (90) does not depend on oj and coinc;i(les with the one for the non-degenerate case. 
Accordingly, the first-order contribution is taken into account just as the solution of the secular equation, i.e. following 
usual techniques developed in quantum mechanics for treating degenerate levels. The influence of quasi-degeneracy 
due to the second-order matrix element V'^^^ is neglected. 

Approximation 3: We calculate the matrix elements of V within the framework of RMBPT. Compared with the 
full ab initio QED calculation the following contributions are missing: 1) negative-energy intermediate states, 2) 
crossed-photon interaction, 3) rigorous treatment of retardation effects. As mentioned above within the framework 
of RMBPT no contribution due to reference states (for two-photon exchange) arises. Accordingly, the energies of the 
levels just follow as solutions of the secular equation. 

Approximation 4: Only the matrix elements of F'-^-' are evaluated within the framework of RMBPT. According to 
the comment made on approximation 2, this also follows the quantum mechanical treatment for quasi-degeneracy. 

Approximation 5: We neglect the imaginary part of the elements of the matrix V. The matrix V defined in Eq. (73) 
is a complex one. Although the energy of the level is defined as the real part of the diagonalized matrix V{w) at the 
point of the resonance {oj = u''''^), the imaginary part of the elements of the matrix V (non-diagonal) contributes to 
the energy. 

The results in Table 1 demonstrate that for the 2^Pi, 2^Pi levels a complete ab initio QED theory for describing 
the quasi-degeneracy has to be employed only when going beyond the level of second-order corrections. For Z < 30 
approximation 3 provides an accuracy of about 1% at the level of second-order perturbation theory. Accordingly, the 
inaccuracy can be referred to corrections of third order. For Z < 60 ~ 70 approximation 4 leads to an inaccuracy 
comparable in magnitude with the corrections of third order. For Z > 60 ~ 70 the effect of quasi-degeneracy decreases 
definitely to the level of third-order corrections. Consequently, approximation 2 can be employed for high-.Z systems. 
For Z > 80 the quasi-degeneracy becomes completely negligible, i.e., it will be sufficient to employ approximation 
1. The contribution of imaginary parts of the matrix elements V to the energy levels appears as an effect of quasi- 
degeneracy, which originates completely from QED. It is perceptible only for high Z > 70, which also reveals that the 
neglect of the imaginary part of y^^' has been legitimate. 

In Table II we present the data for the total energies of the 2^Pi and 2'^Pi two-electron configurations, respectively. 
The numbers present the ionization energy of the 2p-electron with the opposite sign. These data are compared with 
the results obtained by Plante et al. [2] and Drake [1]. Two different approximate methods have been employed in 
these works: the "relativistic all-order theory" (AO) [2] and "the unified theory" [1]. The latter methods accoimt 
approximately for QED effects, such as retardation, crossed-photon graphs and negative-energy intermediate states, 
while taking into account partially higher-orders of the perturbation theory. In the present work the photon exchange 
is taken into account up to the second order. The self-energy (SE) and vacuum-polarization (VP) corrections are 
included only in first order. Quantitative results for SE and VP corrections are taken from Refs. [40-42,3]. The SE 
and VP screening corrections, the radiative corrections of the second order and all the corrections of the third and 
higher orders are omitted. We note, that the VP screening corrections for the states considered have been evaluated 
by Artemyev et al. [15], while results for the SE screening corrections are not yet available. Since the SE and VP 
screening corrections partially cancel each other, we do not include the results of [15] in Table II. In Table III we 
present various theoretical and experimental data for 2'^Pi — 2^Pi transition energies. We conclude that the discrepancy 
between our data and those from other results arising for small values of Z is caused by third- and higher-orders of 
the perturbation theory which have not been accounted for in the present paper. For high Z the major inaccuracy is 
due to missing self-energy, vacuum-polarization screening corrections and one-electron radiative corrections of second 
order. 
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TABLE I. Matrix elements of the operator V for the two-electron configurations 2^Pi and 2^Pi (eV). The individual contribu- 
tions for the Dirac-binding energies of 2p-electron states (y*-"-*), the one-photon exchange contribution {V^^"^) and the two-photon 
contribution (V'^^) are compiled, respectively. i?x(2^Pi), i5x(2''Pi) arc the energies of the corresponding configurations, where 
only the photon-exchange contributions are taken into account (neglecting radiative corrections). 



Contribution 




Z = 10 


18 


26 


30 


40 


(ls2pi/2), (ls2pi/2) 
V^°^: (ls2p3/2),(ls2p3/2) 

(ls2pi/2),(ls2p3/2) 


-340.7099 

-340.2556 




-1108.0574 
-1103.2520 



-2325.7285 
-2304.5586 



-3108.3193 
-3070.5057 



-5594.0369 
-5471.5704 



(ls2pi/2 

(ls2j53/2 
(ls2j3i/2 


),(ls2pi/2) 

),(ls2j33/2) 
), (ls2j33/2) 


64.7130 + 0.0007i 
67.6938 - 0.0007i 
4.3418 - 0.0019Z 


117.2696 -K0.0072i 
122.1610 - 0.0072i 
7.6653 - 0.0204i 


171.2105 0.0316i 
177.1674 - 0.0312i 
10.7333 - 0.0887i 


198.9154 + 0.0560i 
204.9502 - 0.0551i 
12.1369 - 0.1571i 


271.1021 +0.1778i 
275.4795 - 0.1724i 
15.1577 - 0.4954i 


(ls2pi/2), (ls2pi/2) 
(ls2p3/2),(ls2p3/2) 
(ls2pi/2), (ls2j93/2) 


-2.7692 
-3.5256 
—1.0727 


-2.8168 
-3.5603 
— 1.0618 


-2.8938 
-3.6142 
— 1.0450 


-2.9439 
-3.6506 
—1.0350 


-3.1082 
-3.7641 
— 1.0008 


£x(2^Pi) 
Sx(23Pi) 




-273.8939 
—280 9596 


-981.1501 
—997 1 059 


-2127.8323 
—2160 5849 


-2866.5171 
—2915 0368 


-5198.2878 
—5327 6102 


A£;x(2^Pi) 
ASx(2^Pi) 


Appr. 1 
Appr. 1 


2.1934 


3.5012 

— o.OUiz 


3.1731 

— O.i / oi 


2.6890 

— z.DoyU 


1.5672 

— l.OD / Z 


A£;x(2^Pi) 
A£x(23pi) 


Appr. 2 
Appr. 2 


-0.7581 
0.7581 


-0.7203 
0.7203 


-0.5599 
0.5599 


-0.4460 
0.4460 


-0.2175 
0.2175 


A£;x(2^Pi) 
A£;x(2-*Pi) 


Appr. 3 
Appr. 3 


0.0001 
0.0000 


0.0019 
-0.0009 


0.0102 
-0.0036 


0.0168 
-0.0056 


0.0571 
-0.0121 


A£x(2^Pi) 
ASx(2=^Pi) 


Appr. 4 
Appr. 4 


0.0000 
0.0000 


0.0002 
0.0000 


0.0009 
0.0009 


-0.0005 
0.0018 


-0.0013 
0.0038 


ASx(2^Pi) 

ASx(23Pi) 


Appr. 5 

Appr. 5 


0.0000 
0.0000 


0.0000 
0.0000 


-0.0001 
0.0001 


-0.0003 
0.0003 


-0.0015 
0.0015 


Sx(2^Pi) - 
This work 
Lindgren et a 


£;x(2"Pi) 
/. [28] 


7.0657 
7.0657 


15.9557 
15.9554 








TABLE I. {Continued.) 


Contribution 




Z = 50 


60 


70 


80 


92 


(ls2pi/2 

(ls2p3/2 
(ls2pi/2 


),(ls2pi/2) 

),(1.52p3/2) 
), (ls2p3/2) 


-8884.368 
-8575.514 



-13062.966 
-12395.463 



-18250.182 
-16948.025 



-24621.409 
-22253.673 



-34211.065 
-29649.834 



(ls2pi/2 

(ls2p3/2 

(ls2pi/2 


),(ls2pi/2) 
),(ls2p3/2) 
), (ls2p3/2) 


348.917 + 0.437j 
347.960 - 0.4152 
17.320 - 1.206i 


434.638 + 0.912i 
422.949 - 0.846i 
18.430 - 2.491i 


531.400 + 1.702i 
501.042 - 1.531't 
18.311 - 4.590i 


643.793 + 2.930i 
582.883 - 2.533i 
16.795 - 7.777j 


809.699 + 5. 178i 
686.996 - 4.210i 
12.929 - 13.440i 


(ls2pi/2 

(ls2p3/2 

(ls2pi/2 


),(l^.-2pi/2) 

), (ls2p3/2) 
), (ls2p3/2) 


-3.333 
-3.915 

-0.955 


-3.635 
-4.105 

-0.893 


-4.038 
-4.339 
-0.801 


-4.585 
-4.628 
-0.771 


-5.531 
-5.053 

-0.683 


Sx(2^Pi) 
Sx(23Pi) 




-8230.604 
-8539.649 


-11976.160 
-12632.423 


-16451.098 
-17723.044 


-21675.333 
-23982.287 


-28967.898 
-33406.890 


A£x(2^Pi) 
A£;x(2='Pi) 


Appr. 1 
Appr. 1 


0.865 
-0.865 


0.460 
-0.460 


0.225 
-0.225 


0.085 
-0.085 


-0.007 
0.007 


ASx(2^Pi) 

ASx(23Pi) 


Appr. 2 

Appr. 2 


-0.102 
0.102 


-0.049 
0.049 


-0.023 
0.023 


-0.011 
0.011 


-0.004 
0.004 


A£x(2^Pi) 
A£;x(23Pi) 


Appr. 3 
Appr. 3 


0.153 
-0.014 


0.358 
-0.007 


0.759 
0.021 


1.461 
0.096 


2.931 
0.323 


A£;x(2^Pi) 
ASx(23Pi) 


Appr. 4 
Appr. 4 


-0.004 
0.011 


-0.009 
0.025 


-0.015 
0.050 


-0.035 
0.094 


-0.075 
0.181 


A£x(2^Pi) 
ASx(23Pi) 


Appr. 5 
Appr. 5 


-0.004 
0.004 


-0.009 
0.009 


-0.016 
0.016 


-0.026 
0.026 


-0.041 
0.041 
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TABLE II. Data for the energies (in eV) of the configurations 2^Pi and Photon-exchange corrections are taken into 

account up to second order in a. Self-energy (SE) and vacuum-polarization (VP) corrections are taken into account only in 
first order. The onc-clectron radiative corrections of order , the SE and VP screening corrections and all the corrections of 
the third and higher orders are omitted. The data are compared with the results of Plante et al. [2] and Drake [1] . 



Contribution 


Z ^ 10 


18 


26 


30 


40 


E{2'Pi), this work 
E{2^Pi), this work 


-273.8936 
-280.9596 


-981.1462 
-997.1079 


-2127.8070 
-2160.5953 


-2866.4666 
-2915.0545 


-5198.0971 
-5327.6408 


E{2'Pi), Plante et al. 
£;(2^Pi), Plante et al. 


-273.8155 

-281.0132 


-981.0966 

-997.1451 


-2127.7712 

-2160.6313 


-2866.4354 

-2915.0924 


-5198.0801 

-5327.6917 


E{2^Pi), Drake 
E{2^Pi), Drake 


-273.8077 
-281.0054 


-981.0832 
-997.1303 


-2127.7515 
-2160.6038 


-2866.4129 
-2915.0556 


-5198.0515 
-5327.6222 



TABLE II. (Continued.) 



Contribution 


Z = 50 


60 


70 


80 


92 


E{2'Pi), this work 


-8230.079 


-11974.951 


-16448.629 


-21670.722 


-28958.974 


E{2^Pi), this work 


-8539.626 


-12632.120 


-17721.892 


-23979.001 


-33397.135 


E{2'Pi), Plante et al. 


-8230.078 


-11974.972 


-16448.682 


-21670.812 


-28959.135 


E{2^Pi), Plante et al. 


-8539.719 


-12632.321 


-17722.363 


-23980.133 


-33400.643 



E{2^Pi), Drake 



-8539.592 



-12632.088 



-17721.944 



-23979.373 



-33398.993 



TABLE III. Theoretical and experimental data for 2'^Pi — 2^Pi transition energies (in eV). 



z 


This work 


Lindgren et al. [28] 


Plante et al. [2] 


Drake [1] 


Experiment 


Ref. 


18 


15.9617 


16.0550 


16.0485 


16.0471 


16.031(0.074) 


[43] 












16.00(0.49) 


[44] 


26 


32.7883 




32.8601 


32.8523 


33.4(0.5) 


[45] 












33.23(0.45) 


[46,45] 


92 


4438.161 




4441.508 


4440.049 


4455.(87.) 


[47] 
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n 



A 



''XA/^ k', e' 



''X/^ k,e 
A 

FIG. 1. The lowest-order amplitude of the photon scattering on the atomic electron within the resonance approximation. The 
double solid line corresponds to bound electrons in the field of the nucleus. The wavy lines with arrows denote the absorption 
or the emission of a photon with momentum k and polarization e. 



n2 = a 



■XA/^ k', e' 



n\ = a 



'XyQ^ k,e 



A 

FIG. 2. First-order self-energy insertion into the photon-scattering graph within the resonance approximation. 



a 



n 




a 

FIG. 3. The second-order electron self-energy correction (so called SESE loop-after-loop) that gives rise to the correction in 
Eq. (25). 
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"3 = "' 



^3?/^ k',e' 



112 7^ 0- 



ni — a 



'\yi^ k ,e 



FIG. 4. Fcynman graph representing the higher order electron self-energy correction within the line profile approach (SESE, 
loop-after-loop, irreducible). 



k',e' 



k,e 

a 

FIG. 5. Lowest-order amplitudes for photon scattering on a two-electron ion in its ground state A within the resonance 
approximation. The ground state A is represented in terms of non-interacting Dirac electrons. 



k',e' 



ni 



k , e 



n2 



A 



FIG. 6. The lowest-order amplitude for the photon scattering on a two-electron ion in the ground state A- In the ground 
state the interelectron interaction is taken into account. 



FIG. 7. Feynman graph, describing the first-order interelectron interaction. The double solid lines correspond to bound 
electrons in the field of the nucleus, the wavy line corresponds to the exchange of virtual Coulomb and Breit (transverse) 
photons. For a = a and b' = b the graph is called "direct" , and for a = b, b' = a it is called "exchange" graph, respectively. 
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k',e' 
na = a' 



n\ = a 
k,e 



Ui = b' 



n2 = b 



FIG. 8. First-order of interelectron interaction correction to the amplitude of the photon scattering on a two-electron ion 
within the resonance approximation. 



b' a' 



n2 ni 




(a) 

FIG. 9. Fcynman graphs describing the second-order interelectron interaction. The graph (a) is called "box" and the graph 
(b) is called "cross". Notations are the same as in Fig. 7. The summation over intermediate states is indicated by ni,n2. 





(a) (b) 
FIG. 10. Second-order of interelectron interaction correction to the amplitude of the photon scattering on a two-electron ion 
within the resonance approximation. Graph (a) represents the contribution of the "box" graph to the scattering amplitude, 
and graph (b) the contribution of the "cross" graph, respectively . 



a' b' 



"3 

ni 



n4 

n2 



a b 

FIG. 11. The third-order "box" Feynman graph. The notations are the same as in Fig. 7. Here the wavy lines with the cross 
represent the sum of the Coulomb and unretarded Breit interaction. 
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